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ABSTRACT 

Previous observations of quasar host halos at z ~ 2 have uncovered large quantities 
of cool gas that exceed what is found around inactive galaxies of both lower and 
higher masses. To better understand the source of this excess cool gas, we compiled an 
exhaustive sample of 195 quasars at z « 1 with constraints on chemically enriched, cool 
gas traced by Mg II absorption in background quasar spectra from the Sloan Digital 
Sky Survey. This quasar sample spans a broad range of luminosities from Lboi = 
1044.4 io46-8ergs“^ and allows an investigation of whether halo gas properties are 
connected with quasar properties. We find a strong correlation between luminosity and 
cool gas covering fraction. In particular, low-luminosity quasars exhibit a mean gas 
covering fraction comparable to inactive galaxies of similar masses, but more luminous 
quasars exhibit excess cool gas approaching what is reported previously at z ~ 2. 
Moreover, 30 — 40% of the Mg II absorption occurs at radial velocities of |Ar;| > 300 
km s~^ from the quasar, inconsistent with gas bound to a typical quasar host halo. The 
large velocity offsets and observed luminosity dependence of the cool gas near quasars 
can be explained if the gas arises from: (I) neighboring halos correlated through large- 
scale structure at Mpc scales, (2) feedback from luminous quasars, or (3) debris from 
the mergers thought to trigger luminous quasars. The first of these scenarios is in 
tension with the lack of correlation between quasar luminosity and clustering while the 
latter two make distinct predictions that can be tested with additional observations. 
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1 INTRODUCTION 

In order to reproduce the low stellar-to-halo mass ra¬ 
tios of high mass galaxies (e.g. Conroy & Wechsler 2009; 
Behroozi et al. 2013; Kravtsov et al. 2014) both semi- 
analytic (e.g. Benson et al. 2003) and hydrodynamic sim¬ 
ulations of galaxy evolution must incorporate strong feed¬ 
back from active galactic nuclei (AGN) and quasars (for a 
recent review, see Kravtsov & Borgani 2012). Moreover, hy¬ 
drodynamic simulations of galaxy evolution that incorporate 
AGN feedback are better able to reproduce the properties 
of galaxies observed in emission (e.g. Springel et al. 2005; Di 
Matteo et al. 2005; Sijacki et al. 2007; Schaye et al. 2010; Ga¬ 
bor & Dave 2012; Vogelsberger et al. 2013; Li & Bryan 2014). 
While promising, these implementations of quasar feedback 
are subject to significant systematic uncertainties driven by 
poorly constrained aspects of quasar and AGN physics, and 
direct observations of quasar feedback are available for only 
a small number of systems (for a recent review, see Fabian 
2012 ). 

The low-density gas of the circum-galactic medium 
* E-mail: seanjohnson@uchicago.edu 


(CGM) provides a sensitive laboratory for discriminating 
between possible feedback models, and different feedback 
prescriptions result in order-of-magnitude changes in CGM 
observables predicted in hydrodynamic simulations of galax¬ 
ies (Hummels et al. 2013; Ford et al. 2013; Shen et al. 2013a; 
Agertz & Kravtsov 2015; Suresh et al. 2015). The density of 
the CGM is nearly always too low to be studied in emission 
with existing facilities, but significant progress can be made 
by studying the gas in absorption when bright background 
objects are serendipitously found at low projected distances, 
d, from foreground galaxies. 

Over the last decade, large samples totaling nearly one 
thousand galaxies at 2 « 0 to 2 with constraints on extended 
gas properties from absorption spectroscopy have been as¬ 
sembled using a combination of ground and space-based tele¬ 
scopes. In particular, observations of the HI Lyman series 
(e.g. Chen et al. 1998; Tripp et al. 1998; Wakker & Savage 
2009; Stocke et al. 2013; Rudie et al. 2013; Tumlinson et al. 
2013) and doublet transitions due to heavy element ions such 
as Mgll (e.g. Bowen et al. 1995; Chen et al. 2010a; Gauthier 
et al. 2010; Bordoloi et al. 2011; Lovegrove & Simcoe 2011), 
CIV (e.g. Chen et al. 2001; Borthakur et al. 2013; Liang 
& Chen 2014; Bordoloi et al. 2014), and OVI (e.g. Chen 
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& Mulchaey 2009; Wakker & Savage 2009; Prochaska et al. 
2011; Tumlinson et al. 2011; Mathes et al. 2014; Stocke et al. 
2014; Turner et al. 2014; Johnson et al. 2015) have been par¬ 
ticularly fruitful. 

Recently, observing campaigns leveraging large imag¬ 
ing and spectroscopic surveys such as the Sloan Digital Sky 
Survey (SDSS; York et al. 2000) and dedicated campaigns 
on large, ground-based telescopes have extended the study 
of the relationship between galaxy and halo gas properties 
to galaxies hosting quasars using projected quasar-quasar 
pairs. Foreground quasars at both z ~ I and z > 2 exhibit a 
high incidence of optically thick, metal-enriched absorption 
systems traced by HI Lya, CII, CIV, and MgII absorption 
along the transverse direction at d < 300 kpc but low in¬ 
cidence along the foreground quasar sightline itself (Bowen 
et al. 2006; Hennawi et al. 2006; Hennawi & Prochaska 2007; 
Farina et al. 2013; Hennawi & Prochaska 2013; Prochaska 
et al. 2013; Farina et al. 2014; Prochaska et al. 2014). This 
contrast indicates that the ionizing emission from quasars is 
highly anisotropic, in qualitative agreement with the unified 
theory of AGN (e.g. Antonucci 1993; Netzer 2015). 

The high incidence of optically thick gas at d < 300 
kpc from quasars at z « 2 is in significant excess relative to 
that found for inactive galaxies both at z < 1 (Chen et al. 
2010a; Gauthier et al. 2010; Lovegrove & Simcoe 2011) and 
at z fz 2 (Rudie et al. 2012). Insights into the possible ori¬ 
gin of this excess gas can be gained from state-of-the-art 
hydrodynamic simulations that include stellar but not AGN 
feedback (e.g. Fumagalli et al. 2014). In particular, simu¬ 
lations with stellar feedback are able to comfortably repro¬ 
duce the observed incidence of Lyman-limit systems around 
inactive, Lyman-break galaxies at z ~ 2 but under-predict 
the incidence around quasar hosts by more than a factor of 
two (Faucher-Giguere et al. 2015). This discrepancy suggests 
that AGN feedback may be responsible for the excess ob¬ 
served in absorption around quasar hosts. Alternatively, the 
discrepancy could be the result of an inability of the simula¬ 
tions to resolve the physical scales relevant to the formation 
or survival of cool gas clouds in massive halos (Fumagalli 
et al. 2014; Meiksin et al. 2015) or if a substantial portion 
of the gas arises in less luminous galaxies neighboring the 
quasar hosts (Rahmati et al. 2015; Suresh et al. 2015). 

To better understand the relationship between AGN ac¬ 
tivity and halo gas, we searched the SDSS Data Release 
12 (Eisenstein et al. 2011; Alam et al. 2015) and compiled 
an exhaustive sample of 195 quasars at z ~ 1 with con¬ 
straints on Mg H absorption at d < 300 kpc from background 
quasars. This large dataset enables a search for correlations 
between quasar properties and extended, cool circumgalac- 
tic gas around the quasar hosts. 

The paper proceeds as follows: In Section 2 we describe 
the quasar sample and corresponding absorption-line mea¬ 
surements. In Section 3, we characterize the absorption as a 
function of projected distance, redshift, and quasar luminos¬ 
ity. In Section 4, we discuss the implications of our findings. 

Throughout the paper, we adopt a A-cosmology with 
Hm = 0.3, Oa = 0.7, and Ho = 70kms“^ Mpc“^. All mag¬ 
nitudes are in the AB system (Oke & Gunn 1983) and 
corrected for foreground Milky Way extinction following 
Schlafly & Finkbeiner (2011). 


2 QUASAR SAMPLE 

To compile a large sample of foreground-background 
quasar^ pairs with constraints on Mgll absorption, we 
retrieved a list of the 395,281 quasars classified by the 
SDSS-HI automated classification and redshift measurement 
pipeline (Bolton et al. 2012) as of Data Release 12. From 
this sample, we selected foreground-background quasar pairs 
that meet the following criteria: 

(i) The projected distance between the foreground and 
background quasars at the foreground quasar redshift satis¬ 
fies d < 300 kpc. 

(ii) The velocity difference between the foreground and 
background quasars satisfies Au( 2 b,«f) < —10,000 kms“^. 

(iii) The expected wavelengths of the Mg H doublet at the 
redshift of the foreground quasar is covered by the SDSS 
spectrum of the background quasar. 

(iv) The expected wavelengths of the Mg H doublet is out¬ 
side of the Lya forest in the background quasar spectrum, 
and 

(v) The signal-to-noise ratio in the background quasar 
spectrum is sufficient to detect a moderate strength Mgll 
absorption system with rest-frame equivalent width of 
114(2796) > 0.3 A at 3-<j significance at the foreground 
quasar redshift (a signal-to-noise ratio of 20 per SDSS reso¬ 
lution element). 

The upper limit on the projected distance was chosen to 
correspond to the expected virial radii of quasar host halos 
at 2 > 0.4. The requirement that Aw < —10,000 was chosen 
to avoid confusion with gas outflowing from the background 
quasar (e.g. Wild et al. 2008). The requirement that the 
wavelengths of the Mgll doublet at the foreground quasar 
redshift fall outside of the Lya forest in the background 
quasar spectrum was chosen to ensure that absorption at¬ 
tributed to Mgll is not the result of coincidental Lya ab¬ 
sorption systems. Finally, the requirement that the signal- 
to-noise ratio in the background quasar spectrum is suffi¬ 
cient to detect absorption systems of 114(2796) > 0.3 A was 
chosen to correspond to the typical sensitivities of existing 
studies of Mgll absorption around quasars at 2 « 1 (e.g. 
Farina et al. 2013, 2014). 

The search yielded a sample of 195 foreground- 
background quasar pairs which we visually inspected to en¬ 
sure that the automated classifications as broad-line (Type 
1) quasars and redshifts from the SDSS database are ro¬ 
bust. For each of the 195 quasar pairs, we then measured 
foreground quasar properties and Mgll absorption proper¬ 
ties as described in Sections 2.1 and 2.2 respectively. The 
sample is summarized in Table 1. 

2.1 Foreground quasar redshifts and luminosities 

Quasar redshifts from the SDSS pipeline are biased by 
~ 600 kms“^ due to asymmetric, blue-shifted emission pro¬ 
files of broad-line region emitting gas (Gaskell 1982; Tytler 
& Fan 1992; Richards et al. 2002; Hewett & Wild 2010). This 

^ Throughout this paper, we refer to objects with broad-line 
emission and power-law dominated continua as quasars irrespec¬ 
tive of the sub-division of AGN into quasars and Seyferts based 
on luminosity. 
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Table 1. Summary of quasar and absorption properties. The full table is available in the on-line version of the paper. 


Name 


Foreground quasar 


Ae 

(") 

d 

(kpc) 

Absorption 

properties 


Right Ascension 
(J2000) 

Declination 

(J2000) 

2=- 

log Tbol/ergs ^ 

VFr(2796)“^ 

(A) 

Au 

(kms“^) 

Ref.'^ 

J0954-I-3734 

09:54:54.70 

+37:34:19.7 

1.544 

46.3 

3.1 

26 

1.10 + 0.17 

660 

1 

J0836-I-4841 

08:36:49.40 

+48:41:50.0 

0.657 

45.7 

4.1 

28 

1.90 + 0.11 

-130 

1 

J0842-1-4733 

08:42:57.37 

+47:33:42.6 

1.556 

46.5 

3.4 

29 

3.70 + 0.30 

-790 

4 

J1550-I-1120 

15:50:43.59 

+11:20:47.5 

0.4358 

46.0 

5.2 

29 

< 0.25 

n/a 

4 

J1106-I-4635 

11:06:17.17 

+46:35:24.5 

1.602 

46.5 

4.4 

37 

1.12 + 0.15 

-130 

4 

J1108-I-3306 

11:08:07.90 

+33:06:11.3 

1.502 

46.5 

5.5 

46 

5.75 + 0.26 

810 

4 

J0938-I-5317 

09:38:04.21 

+53:17:43.9 

2.063 

45.9 

5.7 

47 

0.48 + 0.04 

310 

3 

J2312-I-1444 

23:12:52.70 

+14:44:58.6 

0.7678 

45.2 

6.4 

47 

0.39 + 0.12 

-100 

1 

J1427-0121 

14:27:58.88 

-01:21:30.3 

2.281 

46.6 

6.2 

51 

0.45 + 0.02 

-80 

3 

J0909-I-1629 

09:09:57.08 

+16:29:06.5 

0.7275 

45.2 

7.1 

51 

< 0.15 

n/a 
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Notes 

Quasars with redshifts from narrow [OII] or [O III] are shown to four decimal places and three decimal places otherwise. 
For non-detections, we report 3-a upper limits integrated over a 250 kms“^ velocity interval. 

Reference: 1 —>• Bowen et al. (2006), 2 ^ Farina et al. (2013, 2014), 3 —>■ Prochaska et al. (2014), 4 ^ This work 


bias exceeds the velocity difference typically found between 
galaxies and associated Mg II absorption at d < 150 kpc 
(Chen et al. 2010a) so more accurate redshifts are required. 
When available, we adopted redshifts from Hewett & Wild 
(2010) and otherwise, we calculated the quasar redshifts us¬ 
ing the template and cross-correlation techniques described 
in Hewett & Wild (2010). In addition, for quasars with nar¬ 
row [OII] or [O III] emission, we measured the redshifts 
from these narrow emission lines by fitting gaussian pro¬ 
files and adopting the rest-frame line centroid wavelengths 
from Hewett & Wild (2010). 

Uncertainties in narrow-line based redshifts due to cen¬ 
troid uncertainties are typically small, cr « 30 kms“^, but 
[OHI] redshifts can be blueshifted due to outflows in the 
narrow-line region (Boroson 2005). To evaluate the signifi¬ 
cance of this bias, we compared [O H] and [O HI] emission 
redshifts for quasars within our SDSS sample and found a 
mean bias of Au(2;oii, zoiii) = —30 kms“^ with a l-cr scat¬ 
ter of 70 kms“^. 

Broad-line based cross-correlation redshifts exhibit 
larger scatter with significant, non-Gaussian wings in 
the redshift error distribution due to population variance 
in broad-line profiles. To evaluate the uncertainties in 
our broad-line based redshifts, we remeasured the cross¬ 
correlation redshifts for quasars with narrow-line emis¬ 
sion but with the narrow-lines masked during the cross¬ 
correlation. We then compared these broad-line redshifts 
with the narrow-line redshifts. The bias in the broad-line 
based cross-correlation redshifts is consistent with zero and 
the 68%, 95%, 99%, and 99.7% uncertainties correspond to 
150, 370,1000, and 1500 kms“^ respectively. 

In addition to improved redshifts, we measured quasar 
bolometric luminosities based on monochromatic, contin¬ 
uum luminosity measurements and the bolometric correc¬ 
tions from Richards et al. (2006). The monochromatic lu¬ 
minosities were estimated by fitting a power-law continuum 
plus Fell template (Vestergaard & Wilkes 2001) model to 
line-free continuum regions of the SDSS quasar spectra with 
rest-frame wavelength intervals of Ar = 1350 — 1360,1445 — 
1465,1700 - 1705, 2155 - 2400, 2480 - 2675, and 2925 - 3500 
A. The uncertainties in the quasar bolometric luminosities 


are dominated by sample variance in quasar spectral energy 
distributions which result in errors of « 0.3 dex (see discus¬ 
sion in Richards et al. 2006). 

In order to include previously reported quasars with ab¬ 
sorption line constraints at d < 300 kpc, we measured bolo¬ 
metric luminosities for quasars from Bowen et al. (2006); 
Farina et al. (2013, 2014) and Prochaska et al. (2014) with 
available SDSS spectra using the same techniques. For the 
quasars from Farina et al. (2013) and Farina et al. (2014) 
without SDSS spectra we adopted the bolometric luminosity 
estimates reported by Farina et al. (2013) and Farina et al. 
(2014) which were measured using a similar continuum fit¬ 
ting procedure. For those quasars from Bowen et al. (2006) 
and Prochaska et al. (2014) without public spectra, we es¬ 
timated luminosities by finding close matches in redshift 
and apparent magnitude space in the quasar sample from 
Shen et al. (2011) which has monochromatic luminosities 
estimated from continuum fitting. These broad-band based 
luminosities reproduce spectroscopic measurements with a 
standard deviation of 0.2 dex. 

The redshift, luminosity, and projected distance distri¬ 
butions of the resulting quasar catalogs are displayed in the 
Figure 1. The samples from Bowen et al. (2006), Farina et al. 
(2013), Farina et al. (2014), and this work are character¬ 
ized by mean redshifts of (z) ~ 1 while the sample from 
Prochaska et al. (2014) is characterized by (z) = 2.2. The 
quasar sample presented in this paper spans a luminosity 
range of log Lboi/ergs“^ = 44.4 to 46.8 (see the left panel 
of Figure 1), extending the study of gas in absorption around 
quasars to lower luminosities than were previously available. 

2.2 Absorption-line measurements 

In order to place constraints on the transverse Mg H ab¬ 
sorption near the foreground quasars, we visually searched 
the background quasar spectra for possible Mgll A2796 
absorption within ]Au] < 1500 kms“^ of the foreground 
quasar redshift. The large search window was chosen to en¬ 
sure that uncertainties in broad-line based quasar redshifts 
do not bias the results of this study. In addition, this ve¬ 
locity search window is the same as was used by Prochaska 
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Figure 1. left: Bolometric luminosity of the foreground quasars versus foreground quasar redshift. The horizontal dotted line marks 
log Ljjoi/ergs"^ = 45.5, the division between the luminous and low-luminosity quasars used throughout the paper, right: Projected 
distance (top) and redshift (bottom) histograms. In all three panels, quasars from Bowen et al. (2006), Farina et al. (2013) or Farina et al. 
(2014), Prochaska et al. (2014), and this work are displayed in cyan, magenta, red, and black respectively, with symbols and line-styles 
indicated in the legend. The quasars from Bowen et al. (2006), Farina et al. (2013), Farina et al. (2014), and this work have available 
constraints on the transverse Mg II absorption from the spectra of background quasars. Mg II absorption constraints are unavailable for 
the quasar sample from Prochaska et al. (2014) due to the higher redshift range. For these quasars, CII rather than Mg II serves as a 
signature of cool, high HI column density gas. 


et al. (2014) which simplihes the comparison between that 
work and our z « 1 sample. 

When absorption was identified, we measured the ab¬ 
sorption equivalent width and centroid in the following 
steps. First, we locally fit the continuum by defining feature- 
free continuum regions at lower and higher wavelengths and 
fit an outlier-resistant line to these continuum regions. We 
then continuum normalized the flux and error arrays and 
calculated the rest-frame equivalent width and equivalent 
width error by direct integration over a user-defined inter¬ 
val. We accepted the absorption according to a 3 -(j detec¬ 
tion threshold and additionally required that the identifi¬ 
cation of the absorption feature as Mg II A2796 be con¬ 
firmed by the presence of at least one other absorption fea¬ 
ture at the same redshift (e.g. Mg II A2803, Mg I A2852, 
or Fell AA2344, 2374, 2382, 2586, 2600). Finally, we mea¬ 
sured the absorber centroid and full-width-at-half-maximum 
(FWHM) by fitting a gaussian profile to the Mg II A2796 
absorption and used the centroid to measure the absorber 
redshift. In a few cases, multiple absorption components are 
found at |Au| < 1500 kms“^. For these quasars, we report 
the total Mg II A2796 equivalent width and report Av for 
the strongest component. Each of these cases is discussed 
further in Section 3.2. 

In the case of non-detections, we placed 3 -ct upper limits 
on the Mg II A2796 equivalent width by integrating the con¬ 
tinuum normalized error array over a 250 kms“^ window 
centered at the foreground quasar redshift. The 250 kms“^ 
window corresponds to the mean FWHM of Mg II A2796 de¬ 
tections in our sample. Two example quasars probed in ab¬ 


sorption, one with detected Mg II absorption and one with¬ 
out, are showcased in Figure 2 to demonstrate the data qual¬ 
ity afforded by the SDSS spectra. 

In one case, J1250—0105, we identified possible Mg II 
A2796 absorption with 3-a significance but were unable to 
confirm the identification by detection of another line. For 
this quasar, we placed a limit on the Mg II A2796 equivalent 
width by measuring the limit on Mg II A2803 and multiplied 
by two (the ratio of the oscillator strengths of the two tran¬ 
sitions) . We note that the results presented in this paper do 
not change if the possible Mg II A2796 absorption for this 
quasar is treated as a detection. 

In a few cases, we detected absorption from the Mg II 
doublet but with the A2796 member contaminated by heavy 
element absorption from another absorption system at a dis¬ 
tinct redshift. For these objects, we fitted the MgII and 
contaminating absorption profiles using the VPFIT package 
(Carswell et al. 1987; Carswell & Webb 2014) and measured 
Mg II A2796 equivalent width from the model fit. 

Finally, we added absorption-line measurements from 
the foreground-background quasar samples previously re¬ 
ported in the literature in Bowen et al. (2006); Farina et al. 
(2013, 2014), and Prochaska et al. (2014). The sample from 
Bowen et al. (2006) contains four foreground quasars probed 
at d = 26,29,47, and 98 kpc all of which are detected 
in Mg II absorption with equivalent widths ranging from 
114(2796) = 0.4 to 1.9 A. We included quasars from Fa¬ 
rina et al. (2013) and Farina et al. (2014) after converting 
their 2-a detection limits to 3 -(t, applying a uniform cut re¬ 
quiring sufficient signal in the background quasar spectrum 
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Figure 2. Two example quasars probed in Mg II absorption. The top panels display a quasar detected in Mg II absorption at d = 72 kpc 
and the bottom panels display a quasar without detected Mg II absorption at d = 52 kpc. The left panels display the SDSS spectra of 
the foreground quasars. The middle panels display the SDSS spectra of the background quasars with the expected positions of the Mg II 
doublet members at the systemic redshift of the foreground quasar marked by vertical dashed lines. In both the left and middle panels, 
the flux and error arrays are shown in black histogram and dotted red line respectively. The zero-flux level is marked by a horizontal, 
black dashed line. The SDSS image of the fields around the quasars are shown on the right with the position of the foreground and 
background quasars labelled as Qf and Qb respectively. The purple line at the top left of each image is b" in length. 


to detect an absorption system of W7(2796) > 0.3 A, and 
restricting the MgII search window to |At;| < 1500 kms”'^. 

The quasar sample from Prochaska et al. (2014) spans 
a redshift range oi z = 1.8 — 3.8, and consequently Mg II 
absorption constraints are not available with the existing 
optical spectra. Mg II absorption traces cool (T ~ 10'* K) 
gas (Bergeron & Stasinska 1986) with high HI column 
densities of ^(HI) « 10*® — 10^^ (Rao et al. 2006). At 
2 ; > 2, CII A1334 absorption can serve as an alternative 
tracer of such cool, high column density gas. Neutral Mag¬ 
nesium and Carbon share similar ionization potential and 
we therefore expect the ratio of Mg II to CII equivalent 
width to approximately follow the rest-frame wavelength ra¬ 
tios, W7(2796)/Wr(1334) = 2.09 for saturated systems with 
line-widths dominated by non-thermal broadening. Empir¬ 
ically, low-redshift galaxies with constraints on both Mg II 
and CII equivalent widths from Werk et al. (2013) exhibit 
a mean ratio of ITr(2796)/W7(1334) = 1.7. In order to com¬ 
pare the cool gas contents around high redshift quasars with 
the sample in this work, we adopted this mean ratio to con¬ 
vert from the CII equivalent widths reported in Prochaska 
et al. (2014) to expected Mg II equivalent widths (consistent 
with the conversion from Prochaska et al. (2014)). To main¬ 
tain a uniform cut in absorption sensitivity, we restricted 
the quasar sample from Prochaska et al. (2014) to those 
with signal-to-noise in the background spectrum sufficient to 
detect a CII A1334 absorption system of 117(1334) > 0.18 


A at 3 -ct significance. This sensitivity corresponds to the 
117(2796) > 0.3 A sensitivity requirement adopted for the 
z ~ 1 samples with constraints on Mg II absorption. 


3 MG II ABSORPTION NEAR QUASARS 

With the quasar and absorption-line measurements de¬ 
scribed in Section 2 in hand, we characterize the Mg II ab¬ 
sorption found around quasars as a function of projected 
distance, luminosity, and redshift in Section 3.1. In addi¬ 
tion, we explore the kinematics of the absorption systems in 
Section 3.2 and discuss associated Mg II absorption systems 
found along the sightline to the foreground quasar itself in 
Section 3.3. 


3.1 Dependence on projected distance, 

Inminosity, & redshift 

We characterize Mg II equivalent width as a function 
of projected distance in the left panel of Figure 3 and dis¬ 
play the measurements from Bowen et al. (2006); Farina 
et al. (2013, 2014), and Prochaska et al. (2014) for compari¬ 
son. The Mg II equivalent widths exhibit an anti-correlation 
with projected distance that is driven by a decrease in inci¬ 
dence at large projected distances. To verify the significance 
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Figure 3. left: Rest-frame Mg II A2796 equivalent width versus projected distance. Symbols and coloring are as in Figure 1. Non-detections 
are shown as 3-cr upper limits marked by downward arrows and lighter shading, right: Covering fraction for absorption systems with 
Wr(2796) > 0.3 A as a function of projected distance for the quasars probed in Mg II with luminous quasars (log Lboi/ergs~^ > 45.5) 
in black solid line and low-luminosity quasars (log Lboi/ergs“^ < 45.5) in blue dashed line. For comparison, the covering fraction for 
luminous quasars at 2 ps 2.2 based on sample probed in CII from Prochaska et al. (2014) is shown in red dotted line. Finally, we also 
show the expected covering fraction for galaxies with stellar mass log M^/Mq = 11 in green dash-dot line. The vertical error bars mark 
68% confidence intervals calculated from binomial statistics with a fiat prior [k G [0,1]) and the horizontal error bars mark the full range 
of projected distances contributing to each bin. The bins are chosen to span the range of d = 0 — 300 kpc in 100 kpc wide intervals. We 
note that the expected incidence of coincidental Mg II absorption systems unrelated to the quasars within the velocity search window 
is insignificant (/^rand = 0.01) based on the Mg II dN/dzdWx measurement from Zhu ^ Menard (2013). The SDSS quasars are drawn 
from a flux-limited survey and consequently, the luminous quasars occur at a higher median redshift of (2)med = 0.98 compared to 
(2}med = 0.75 for the low-luminosity quasars. The higher gas covering fractions observed for the luminous quasars could therefore be the 
result of luminosity-dependence, redshift evolution, or a combination of the two. The observed correlation between covering fraction and 
luminosity is isolated from the possible effects of redshift evolution in Figure 4. 
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of this anti-correlation between equivalent width and pro¬ 
jected distance, we perform a generalized Kendall-r rank 
correlation test including the non-detections as upper limits 
(Isobe et al. 1986). The Kendall test finds an anti-correlation 
characterized by a coefficient of —0.3 with a high signifi¬ 
cance of p ^ 1%^. The anti-correlation is driven primar¬ 
ily by a decrease in the incidence of absorption systems at 
d > 100 kpc rather than by a decrease in the mean equiv¬ 
alent width of detections. Interestingly, strong absorption 
systems of 144(2796) > 1.0 A are observed at d > 100 kpc, 
in stark contrast to L* galaxies at z < 1.5 (Chen et al. 2010a; 
Lovegrove & Simcoe 2011). 

Using a sample of quasars probed in absorption by back¬ 
ground quasar spectra placing constraints on C If absorp¬ 
tion, Prochaska et al. (2014) found excess absorption for 
quasars at 2 « 2.2 compared to quasars probed in Mg II at 
z ~ 1 (Farina et al. 2013, 2014). Prochaska et al. (2014) in¬ 
terpreted this excess as redshift evolution in the properties of 


^ We note that the measurement uncertainties will tend to bias 
the correlation coefficient toward zero and perform this test to 
estimate the statistical significance of the anti-correlation within 
our sample rather than to measure the coefficient. 


quasar host halos. However, the quasars in the 2 « 2.2 sam¬ 
ple from Prochaska et al. (2014) are characterized by high lu¬ 
minosities (see left panel of Figure 1), and it is not clear that 
they form a representative sample for the entire quasar pop¬ 
ulation. To investigate this possibility, we combine our SDSS 
quasars with those reported in the literature with Mg II con¬ 
straints (Bowen et al. 2006; Farina et al. 2013, 2014) and 
divide the resulting combined Mg II sample at the median 
luminosity into luminous quasars of log Lboi/ergs”'^ >45.5 
and low-luminosity ones of log Lboi/ergs“^ < 45.5. In do¬ 
ing so, we restrict the combined Mg II quasar sample to 
those with 2 < 1.8 where the sample follows the well defined 
redshift-luminosity trend seen in Figure 1. We note that this 
luminosity division also corresponds to the traditional divi¬ 
sion of AGN into quasars and Seyfert galaxies (e.g. Schmidt 
& Green 1983). 

Both the luminous and low-luminosity quasars exhibit 
decreasing covering fraction as a function of projected dis¬ 
tance as shown in the right panel of Figure 3. In particular, 
the low-luminosity quasar covering fraction decreases from 
KMgii = 0.50±0.14 at d < 100 kpc to ACMgii = O.lOlg;®® and 
KMgii = 0.021 q'qj at d = 100 — 200 and 200 — 300 kpc re¬ 
spectively. The Mg II covering fractions for low-luminosity 
quasars are consistent within uncertainties with expecta- 
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tions for inactive galaxies with stellar masses log M^/Mq « 
11 based on the covering fraction and stellar mass scaling 
relation measurements from Chen et al. (2010a) and Chen 
et al. (2010b)^. A stellar mass of log M^/Mq « 11 is cho¬ 
sen for this comparison because it is the expected stellar 
mass of galaxies with halo masses similar to quasar hosts 
(Shen et al. 2013b) based on the stellar-to-halo mass rela¬ 
tion from Behroozi et al. (2013). On the other hand, the 
luminous quasars exhibit enhanced covering fraction that 
are a factor of « 2 or more higher than those of the low- 
luminosity quasars over the full range of projected dis¬ 
tances studied here with fvMgii = O.SSlg'jg at d < 100 kpc, 
KMg II = 0.42 ± 0.07 at d = 100 — 200, and Kiug ii = O.ldjllQ'gg 
at d = 200 — 300. 

To verify the significance of the excess Mg II absorp¬ 
tion around luminous quasars relative to low-luminosity ones 
over the full range of projected distances studied here, we 
also calculate the covering fraction at d < 300 kpc for 
the low-luminosity quasars and luminous quasars and find 
KMgii = O.IOIq'qj and 0.347Q g 4 respectively. In addition, we 
perform a logrank test (Feigelson & Nelson 1985) comparing 
the distributions of the Mg II equivalent widths for the low- 
luminosity quasars and redshift weighted luminous quasars 
at d < 300 kpc. The logrank tests yields Z/„ = 15.2 ± 3.7, 
which confirms that the luminous quasars exhibit excess 
Mg II absorption relative to low-luminosity quasars with a 
signihcance of p ^ 1%. 

While the low- and high-luminosity quasars exhibit dis¬ 
tinct gas covering fractions, it is possible that the differ¬ 
ence in gas covering fraction is the result of redshift evolu¬ 
tion in quasar host halo gas properties. In particular, the 
SDSS quasar sample is drawn from a flux-limited survey, 
and consequently the luminous quasars in our sample oc¬ 
cur a higher median redshift of (2:)med = 0.98 compared 
to (z)med = 0.75 for the low-luminosity quasars as shown 
in Figure 1. The excess absorption observed around the lu¬ 
minous quasars could therefore be the result of luminosity 
dependence, redshift evolution, or a combination of the two. 
To isolate the luminosity dependence, we restrict the lumi¬ 
nous quasar sample to 2 < 1 in order to match the median 
redshift of the luminous quasar sample to that of the low- 
luminosity quasars. The covering fraction for the luminous 
quasars at 2 < 1 and d < 300 kpc is Kivigii = 0.32to o 5 com¬ 
pared to KMgii = O.IOIq qj for the low-luminosity quasars. A 
logrank test comparing the Mg II equivalent width distribu¬ 
tion of the low-luminosity quasars with that of the luminous 
quasar at 2 < 1 confirms this excess with a significance of 
p 1%. The comparison of the Mg II absorption incidence 
between the luminous and low-luminosity quasars with re¬ 
stricted redshift ranges indicates that the observed correla¬ 
tion between covering fraction and luminosity is not being 
driven by underlying redshift evolution. 

The observed correlation between Mg II absorption cov¬ 
ering fraction and quasar luminosity can explain a significant 
portion of the excess cool gas found for quasars at 2 « 2.2 
relative to those at 2 « 1. Nevertheless, the covering frac- 


^ The galaxy sample from Chen et al. (2010a) is based on a sam¬ 
ple of L Lt galaxies at 2 < 0.2, but the scaling relations in 
Chen et al. (2010b) remain valid at 2 < 2 (Lovegrove & Simcoe 
2011; Chen 2012; Liang & Chen 2014). 


tions for the luminous quasars at 2 « 1 probed in Mg II at 
d = 100 — 200 kpc are somewhat lower than those of the 
2 « 2.2 quasars probed in CII from Prochaska et al. (2014) 
(k — 0.69lg'jg). This possible difference between the z ~ 1 
and 2 « 2 samples could be the result of additional cor¬ 
relation with luminosity, a difference in projected distance 
distribution between the 2 « 1 and 2 « 2.2 samples, or 
evolution in cool gas content of quasar host halos. 

To differentiate these possibilities, we plot the covering 
fraction of Mg II absorption systems with 114(2796) > 0.3 A 
at d < 200 kpc as a function of luminosity for the combined 
Mg II quasar sample after resampling to mimic the flatter 
projected distance distribution of the Prochaska et al. (2014) 
sample in the left panel of Figure 4. The resampling is ac¬ 
complished by randomly resampling the quasars probed in 
Mg II without replacement to construct a maximal possible 
sample with the projected distance distribution mimicking 
that of Prochaska et al. (2014). The random resampling is 
repeated a large number of times to measure the mean cover¬ 
ing fraction. Uncertainties in the mean covering fraction are 
calculated from binomial statistics with the sample size of 
the largest possible subsample with fiat projected distance 
distribution. 

The low-luminosity quasars (log (I/boi/ergs“^) = 45.2) 
exhibit a mean covering fraction of Kivigii = 0.25lg'gg while 
the luminous ones (log (Lboi/erg s“^) = 45.9) exhibit a 
covering fraction of and fCMgii = O.GIq q® at d < 200 
kpc, indicating a strong correlation between quasar lumi¬ 
nosity and the incidence of extended Mg II absorbing gas. 
A similar trend of increasing covering fraction with lumi¬ 
nosity is observed when the sample is restricted to 2 < 1 
(see the blue points in Figure 4). The trend of increas¬ 
ing covering fraction with luminosity found for the 2 « 1 
quasars connects smoothly with the high covering frac¬ 
tions (k = 0.767o;i 3) found for high-luminosity quasars 
(log (Lboi/ergs“^) = 46.2) at 2 « 2.2 from Prochaska et al. 
(2014) (see the left panel of Figure 4). We note that in cal¬ 
culating the covering fraction for quasars from Prochaska 
et al. (2014), we restrict the sample to 2 < 2.5 where the 
Prochaska et al. (2014) sample closely follows the redshift- 
luminosity trend seen in the 2 « 1 samples (see the left panel 
of Figure 1). 

With the correlation between Mg II absorption inci¬ 
dence and quasar luminosity in mind, we search for any 
additional dependence of Mg II absorption on redshift. To 
do so, we plot the covering fraction at d < 200 kpc of the lu¬ 
minous quasars as a function of redshift after resampling to 
mimic the projected distance distribution of the Prochaska 
et al. (2014) quasar sample in the right panel of Figure 4. 
At 2 < 1 and 2 = 1 — 2, the luminous quasars probed in 
Mg II exhibit covering fractions of Kivigii = 0.55lg'jg and 
0.647 o'io at d < 200 kpc, consistent with no evolution be¬ 
tween 2 < 1 and 2 = 1 — 2. The corresponding covering 
fraction for the 2 « 2.2 sample from Prochaska et al. (2014) 
is K = 0.76 ±0.12, consistent with the covering fractions 
found for the 2 < 1 quasars at the level of 1.3 ct. The in¬ 
crease from KMgii = 0.53tQ'42 at 2 < 1 to k = 0.76jlQ''j'3 at 
2 « 2.2 is suggestive, but a larger sample of quasars at 2 « 2 
is required to investigate further. In addition, optical spectra 
of background sightlines toward quasars at 2 < 2.5 can pro¬ 
vide constraints on Mg II absorption which would eliminate 
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Figure 4. left Covering fraction at d < 200 kpc of Mg II absorption systems of 114(2796) > 0.3 A as a function bolometric luminosity 
of the foreground quasars. To isolate the correlation with luminosity from possible redshift evolution, we also show the covering fraction 
as a function of bolometric luminosity with the quasar sample restricted to 2 < 1 in blue. This redshift restriction is chosen so that the 
median redshift of the luminous quasars matches that of the low-luminosity quasars, right. Covering fraction at d < 200 for quasars of 
log Luol/ergs"^ > 45.5 as a function of redshift. In both panels, the vertical error bars mark 68% confidence intervals estimated from 
binomial statistics and the horizontal error bars mark the full range contributing to each bin. In both panels, the MgII quasar sample 
is resampled to mimic the projected distance distribution of the z Ri 2.2 sample form Prochaska et al. (2014) (red). 


the additional systematic uncertainty in converting between 
Mg II and CII absorption equivalent widths. 

3.2 Kinematics 

To characterize the kinematics of the Mg II absorption 
systems around quasars, we plot histograms of the radial 
velocity differences between the quasars and Mg II systems 
in the left panel of Figure 5. The radial velocity histograms 
are characterized by a core that is consistent with the veloc¬ 
ity spread found around galaxies (a « 150 kms”'^; Chen 
et al. 2010a), but 29 out of 67 Mg II components (43%) 
are found at radial velocity differences with |Au| > 300 
kms“^ compared to 2 out of 47 galaxies (4%) in Chen et al. 
(2010a). Ten of 29 quasars with narrow-line redshifts and 
detected Mg II absorption are found to have Mg II absorp¬ 
tion at |Au| = 300 — 1500 kms“^ (34%) indicating that 
the majority of the large velocity differences are not the 
result of the uncertainties in broad-line based quasar red- 
shifts (see the top left panel of Figure 5). Such large velocity 
differences are found among both high- and low-luminosity 
quasars, though the sample of low-luminosity quasars with 
Mg II absorption is small (see the bottom left panel of Figure 
5). 

In addition to large radial velocity differences, six 
quasars in the SDSS sample are found to have multiple, 
resolved Mg II absorption complexes spread over hundreds 
to thousands of kms“^ in the background quasar spectrum. 
The complex kinematics found for these six quasars are vi¬ 
sualized in the right panel of Figure 5. The Mg II doublet 


members are separated by a wavelength difference that cor¬ 
responds to Av = 770 kms“^, creating the potential for con¬ 
fusion between absorber kinematics and the weaker doublet 
member. To avoid this confusion in visualizing the systems 
of multiple Mg II absorption complexes, we ht the Mg II ab¬ 
sorption using the VPFIT package and plot the data after 
dividing by the ht for the A2803 doublet member. We note 
that each absorption “complex” in the SDSS spectra would 
likely be resolved into multiple components in high resolu¬ 
tion spectra. 


3.3 Narrow, associated absorption systems 

A few percent of optically selected quasars exhibit nar¬ 
row, “associated” absorption-line systems (zabs ~ Zqso) 
along the sightline to the quasar itself (e.g. Weymann et al. 
1979; Wild et al. 2008; Vanden Berk et al. 2008). Associ¬ 
ated absorption-line systems could arise in gas associated 
with the quasar, the interstellar medium of the quasar host, 
or the interstellar/circum-galactic medium of neighboring 
galaxies, but the ionizing radiation from the quasar is ex¬ 
pected to photoevaporate Mg II absorbing clouds out to dis¬ 
tances of « 1 Mpc along the quasar sightline (Hennawi & 
Prochaska 2007; Wild et al. 2008). A recent comparison be¬ 
tween quasars with narrow, associated Mg II absorption and 
a control sample of redshift and i-band magnitude matched 
quasars without associated absorption found that quasars 
with associated Mg II absorption exhibit slightly enhanced 
dust extinction {E{B — V)~ 0.03) and a 50 — 70% enhance¬ 
ment in [Oil] emission (Shen & Menard 2012). 
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Figure 5. left: Histograms of the radial velocity differences between the quasar systemic redshifts and MgII component redshifts at 
Av < 1500 kms“^. The top panel shows the full MgII sample in black, quasars with narrow-line redshifts in shaded red, and the fit 
to galaxies from Chen et al. (2010a) in green dash-dot line. Vertical lines mark ±300 kms“^ , the expected virial velocity of a typical 
quasar host halo a,t z = 1.0. The bottom panel shows the radial velocity histogram for the high- and low-luminosity quasars in black and 
blue filled histogram respectively, right: Kinematic structure seen in absorption toward six foreground quasars observed with resolved 
velocity structure in SDSS spectra of the background quasar. The panels are labelled with the foreground quasar name, redshift, projected 
distance, and bolometric luminosity. Quasar redshifts from narrow lines are shown to four decimal places and those from broad lines are 
shown to three. To prevent confusion between kinematic structure and the weaker Mg II doublet member, we model the Mg II doublet 
observed in the SDSS spectra and plot the data divided by the model for the A2803 doublet member. 


A search within our foreground quasar sample identifies 
four luminous and six low-luminosity quasars with associ¬ 
ated (I At)I < 1500 kms“^) narrow MgII absorption systems 
along the foreground quasar sightline. Two of these four lu¬ 
minous quasars are found to have Mg II absorption systems 
detected in the background quasar spectra at d = 150 — 300 
kpc, consistent with expectations from the general high- 
luminosity quasar sample. None of the six low-luminosity 
quasars with associated absorption systems are found to 
have Mg II systems in the background quasar spectra at 
d = 180 — 275 kpc, also consistent with expectations from the 
general, low-luminosity quasar sample. Though the sample 
of quasars with associated Mg II absorption and constraints 
on transverse Mg II absorption from background quasars is 
small, these results suggest that the large-scale halo gas con¬ 
tents of quasars with associated absorption are not radically 
different from those found in the general quasar population. 


4 DISCUSSION 

Using a sample of 195 quasars at 2 : ~ 1 with constraints 
on Mg II absorption from background quasars at projected 
distances of d < 300 kpc from the SDSS, we characterized 
the cool gas contents of quasar host halos as a function 
of projected distance, quasar luminosity, and redshift. Our 
main findings are the following: 

(i) Luminous quasars of log Lboi/ergs”'^ > 45.5 ex¬ 
hibit enhanced Mg II absorption relative to low-luminosity 


quasars and inactive galaxies of similar mass at projected 
distances of d < 300 kpc. 

(ii) The absorbing gas near quasars exhibits complex 
kinematics with 30 — 40% of components found at |Au| = 
300 — 1500 kms“^ from the quasar systemic redshift. and 

(iii) The incidence of cool gas absorption around lumi¬ 
nous quasars does not evolve strongly with redshift between 
2 ~ 1 and 2 « 2 . 2 . 

In this section, we discuss the possible origins of the 
extended Mg II absorbing gas near quasars, correlation with 
luminosity, and complex kinematics. Finally, we briefly con¬ 
sider possible avenues of future research. 

Galaxies of L « L* at 2 < 0.2 exhibit high Mg II cover¬ 
ing fractions out to a gaseous radius that is observed to scale 
with stellar mass according to Rms II oc (Chen et al. 

2010b) and comparisons with CGM absorption at higher 
redshift suggest that the low-redshift scaling relations re¬ 
main valid at 2 < 2 (e.g. Lovegrove & Simcoe 2011; Chen 
2012; Liang & Chen 2014). Since quasars are thought to re¬ 
side in massive galaxies of L > L*, previous studies (e.g. 
Farina et al. 2014; Prochaska et al. 2014) suggested that the 
high incidence of CII and Mg II absorption observed at im¬ 
pact parameters of d = 100—200 kpc from quasars are the re¬ 
sult of mass scaling of the CGM. In this scenario, the absorp¬ 
tion traces the “normal” halo gas of inactive galaxies with 
masses similar to quasar hosts. However, the expected Mg II 
absorption incidence from inactive galaxies with masses sim¬ 
ilar to quasar hosts (log M*/M 0 « 11 based on the cluster- 
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ing measurement from Shen et al. (2013b) and stellar-to-halo 
mass relation from Behroozi et al. (2013)) is significantly 
lower than the covering fraction observed around luminous 
quasars (see the right panel of Figure 3). Moreover, quasar 
luminosity and host halo mass are only weakly (if at all) cor¬ 
related (e.g. Shen et al. 2013b) so the observed correlation 
between quasar luminosity and Mg II absorption cannot be 
explained by mass scaling. 

Finally, the large velocity differences observed between 
the quasar and absorber redshifts (30 — 40% of absorption at 
|Au| = 300 — 1500 kms“^) and complex kinematics seen in 
some sightlines (see Figure 5) are inconsistent with gas grav¬ 
itationally bound to a halo with the mean mass of quasar 
hosts (log Mh/Mq = 12.8; Shen et al. 2013b) which have 
inferred halo virial velocities of ~ 300 kms“^. However, 
clustering measurements constrain the mean mass of quasar 
hosts, not the overall host halo mass distribution (see dis¬ 
cussion in White et al. 2012), and the large fraction of 
absorption at |Aii| = 300 — 1500 kms“^ could be explained 
if « 30 — 40% of quasars reside in halos with larger masses. 
To evaluate this possibility, we measure the radial velocity 
differences observed between quasars with narrow-line red- 
shifts and Mg II absorption systems, correct for line-of-sight 
projection by multiplying by \/3, and find a population aver¬ 
aged velocity dispersion of 1000 kms“^ between the quasars 
and absorbing gas. The halo mass corresponding to a virial 
velocity of 1000 kms“^ is log M^/Mq = 14.4. Such mas¬ 
sive halos are three orders-of-magnitude less abundant than 
those of log Mu/Mq = 12.8 (Tinker et al. 2008) making it 
unlikely that a substantial fraction of quasars reside in such 
massive systems. Direct confirmation of this conclusion will 
require redshift surveys to characterize the environments of 
quasars with constraints on halo gas properties from back¬ 
ground sightlines. 

The high Mg II gas covering fractions, correlation with 
quasar luminosity, and kinematics can be explained if the 
Mg II absorbing gas is the result of: 

(i) the CGM of neighboring galaxies at > 1 Mpc scales 
but at projected distances of d < 300 kpc, 

(ii) feedback from luminous quasars, or 

(iii) debris from the galaxy interactions and mergers 
thought to trigger luminous quasars. 

In the following subsections, we discuss each of these possi¬ 
bilities in turn. 

4.1 Mg II absorption from galaxies neighboring 
the quasar host 

Large scale quasar-quasar (e.g. White et al. 2012) and 
quasar-galaxy (e.g. Shen et al. 2013b) clustering measure¬ 
ments indicate that quasars reside in massive halos that 
trace over-dense regions of the Universe. Consequently, the 
CGM of galaxies in neighboring halos at distances of > 1 
Mpc from the quasar but projected distances of < 300 kpc 
could contribute to the Mg II absorption covering fractions 
and explain the large radial velocity differences often ob¬ 
served between the quasar and absorber redshifts. 

To evaluate the expected covering fraction from neigh¬ 
boring galaxies, we use four UV-bright quasars at z = 
0.3 — 0.6 with highly complete redshift surveys of galax¬ 
ies of Tab < 23 from Ghen & Mulchaey (2009); John¬ 


son et al. (2013); and Johnson et al (in prep). These sur¬ 
veys reveal the presence of galaxies with redshifts within 
|Au| < 1500 kms“^ of the quasar redshift at projected dis¬ 
tances of d < 300 kpc that could potentially contribute to 
the Mg II absorption observed around quasars at similar pro¬ 
jected distances (see left panel of Figure 6). We estimate 
the expected contribution to the covering fraction observed 
around quasars from these neighboring galaxies using the 
covering fraction and luminosity scaling measurements from 
Ghen et al. (2010a). In this calculation, we do not include 
galaxies more luminous than Mb = —22.2 since such lumi¬ 
nous galaxies exhibit reduced incidence of Mg II absorption 
compared to galaxies of L « L* (Gauthier et al. 2010). 

The expected Mg II covering fractions from neighboring 
galaxies can explain a significant portion of the Mg II ab¬ 
sorption observed around quasars at d < 300 kpc (see right 
panel of Figure 6). Moreover, neighboring galaxies are often 
at |Au| « 300 — 1500 kms“^ from the quasar, potentially 
explaining the kinematics observed in Mg II absorption. We 
note, however, that the sample of quasars with available 
deep redshift surveys is small (4 quasars) and, consequently, 
the estimated covering fraction from neighboring galaxies 
suffers from significant and poorly quantified sampling un¬ 
certainties. 

If a substantial portion of the Mg II absorption observed 
around quasars traces the CGM of neighboring galaxies in 
other host halos at Mpc scales, then the correlation between 
quasar luminosity and Mg II absorption would imply a cor¬ 
relation between quasar luminosity and quasar-galaxy clus¬ 
tering. However, studies of quasar-quasar and quasar-galaxy 
clustering found no correlation between quasar luminosity 
and clustering on Mpc scales (e.g. Groom et al. 2005; Myers 
et al. 2007; White et al. 2012; Shen et al. 2013b). In ad¬ 
dition, recent observations of projected quasar-photometric 
galaxy pair counts find no evidence for a correlation with 
quasar luminosity (Padmanabhan et al. 2009; Zhang et al. 
2013; Scott et al. 2015). We note that some studies of quasar- 
photometric galaxy pair counts find excess counts at d < 300 
kpc for luminous quasars (Serber et al. 2006), but this ex¬ 
cess is attributed to galaxies in the quasar host halo and the 
excess is not observed on larger scales. 

The lack of correlation between quasar luminosity and 
clustering are in tension with expectations from the observed 
correlation between Mgll covering fraction and quasar lu¬ 
minosity if a substantial portion of the absorption arises in 
neighboring halos on Mpc scales, in apparent contradiction 
with the covering fraction estimate shown in the right panel 
of Figure 6. This discrepancy can be explained by the prox¬ 
imity effect in which the UV light from the quasar photo- 
evaporates the cool halo gas of galaxies out to Mpc scales. 
Indeed, galaxies at Au < 3000 kms“^ from quasars exhibit 
significantly reduced cool halo gas content compared to the 
general galaxy population (e.g. Pascarelle et al. 2001). 

4.2 Mg II absorption due to quasar feedback 

Radio-loud, lobe-dominated quasars (class FRH; Fa- 
naroff & Riley 1974) are known to drive powerful outflows at 
both high and low redshift (e.g. Nesvadba et al. 2008; Fu & 
Stockton 2009). However, only a small fraction quasars are 
radio-loud and lobe-dominated (e.g. Urry & Padovani 1995; 
Ivezic et al. 2002), so radio-mode feedback alone cannot be 
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Figure 6. left: The galaxy environments near a luminous quasar (3C 263 at 2 : = 0.646) with available deep redshift surveys targeting 
galaxies of tab < 23. Galaxies with secure redshifts and within |An| < 1500 kms“^ of the quasar redshift are labelled with An and 
star-forming galaxies are highlighted with dotted black outline. A dashed black circle of 150 kpc in radius centered on the quasar is shown 
for scale. The image is a stack of two 600 second exposures from the Advanced Camera for Surveys on the Hubble Space Telescope (PI: 
Mulchaey; PID: 13024) taken with the F814W filter with offsets between exposures to fill the gap between the ACS detectors. Cosmic ray 
removal is not possible in regions with coverage in only one exposure and these regions are visible as 5” wide stripes with high cosmic ray 
contamination, right: Covering fraction for absorption systems with 144(2796) > 0.3 A as a function of projected distance for luminous 
quasars (black solid line) compared to expectations from a galaxy of logMt/MQ = 11 (green dash-dot line), and expected incidence due 
to neighbors of the quasar host (blue dashed line). The mean expected incidence due to neighbors is based on four UV-bright quasars 
with available deep galaxy redshift surveys and the light blue band represents the estimated uncertainty in the mean. We note, however, 
that with such a small sample size, uncertainty due to sampling are both large and poorly quantified. 




responsible for the observed correlation between Mg II ab¬ 
sorption and luminosity. 

Recent observations of spatially resolved [OIII] emis¬ 
sion around luminous, radio-quiet, obscured (Type 2) 
quasars revealed the presence of surprisingly spherical 
gaseous nebulae that are both spatially (« 30 kpc) and kine¬ 
matically (« 1000 kms“^) extended (Greene et al. 2012; 
Liu et al. 2013a,b). Similar structures are observed in [OIII] 
emission around luminous, radio-quiet, unobscured (Type 
1) quasars (Liu et al. 2014) like those studied in this work. 
These observations are most naturally explained as fast out¬ 
flows with wide opening angles driven by radio-quiet quasars 
(also see Zakamska & Greene 2014). 

The quasar driven outflows observed in emission on 
scales of « 30 kpc may naturally explain the high incidence 
and complex kinematics of Mg II absorption around quasars 
if the outflows persist to larger distances but with gas den¬ 
sities that are too low to be observed in collisionally excited 
emission lines. An outflow with velocity Uout = 1000 kms“^ 
at 15 kpc from the host galaxy nucleus could reach distances 
of « 100 (200) kpc in « 10® (2 x 10®) years assuming that 
the outflow decelerates due to the gravitational potential of 
the host halo. The timescales required to reach 100 — 200 kpc 
from the host halo are comparable to the high end of quasar 
lifetime estimates (10® years; Martini 2004). Moreover, these 


fast, quasar driven outflows are observed to be effective at 
luminosities of log Lboi/ergs“^ > 45.5 (e.g. Veilleux et al. 
2013; Zakamska & Greene 2014) and could therefore explain 
the observed correlation between extended Mg II absorption 
and quasar luminosity. 

In the outflow scenario, the Mg II absorbing gas repre¬ 
sents cool clumps entrained in a hotter outflowing medium 
(e.g. Gosta et al. 2015). Indeed, observations of outflows 
from both quasars and ultra-luminous infrared galaxies re¬ 
veal copious quantities of highly ionized gas traced by CIV 
and O VI absorption (Arav et al. 2013; Martin et al. 2015). 
However, in the one foreground quasar with published high 
resolution background quasar spectroscopy, cool gas traced 
by CII absorption is found spread over 700 kms“^ but lit¬ 
tle absorption from more highly ionized species is found 
(Prochaska & Hennawi 2009). The low ionization state of 
the gas in this system (AI(C II)/A^(C IV) > 10) is inconsis¬ 
tent with expectations from cool gas entrained in a quasar 
driven outflow, suggesting that quasar driven outflows alone 
cannot universally explain the high absorption covering frac¬ 
tions and complex kinematics observed around quasars. 

Alternatively, the gas kinematics and ionization state 
can be understood if feedback from luminous quasars drives 
the pre-existing halo gas to large velocities, far in excess of 
the sound speed. Such fast bulk motions could lead to an en- 
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hanced fraction of gas with sufficiently high densities to cool 
and produce Mg II absorbing clumps. In this way, feedback 
from the quasar could convert halo gas generally observed 
in high ionization states at projected distances less than 
the virial radius (e.g. O VI; Chen & Mulchaey 2009; Tum- 
linson et al. 2011; Johnson et al. 2015) to lower-ionization 
states observable in Mg II absorption. Feedback from lumi¬ 
nous quasars therefore represents a viable scenario for ex¬ 
plaining the cool halo gas properties around quasar hosts 
including kinematics, correlation with luminosity, and en¬ 
hanced covering fractions relative to inactive galaxies of sim¬ 
ilar masses. 


4.3 Mg II absorption from tidal debris 

Luminous quasars are thought to be fueled by gas sup¬ 
plied during galaxy mergers while the gas required to fuel 
less luminous AGN can be supplied by more secular pro¬ 
cesses (e.g. Hopkins & Hernquist 2009). The high incidence 
of absorbing gas around quasars and correlation with lumi¬ 
nosity can be explained if the gas traces debris from the 
galaxy mergers that can trigger luminous AGN. Tidal de¬ 
bris can obtain velocities that exceed the escape velocity 
of the host halo (e.g. Toomre & Toomre 1972), possibly 
explaining the kinematics observed in absorption. The low 
ionization state (A^(C II)/A1(G IV) > 10) observed in the 
one available quasar with high resolution background spec¬ 
troscopy (Prochaska & Hennawi 2009) is consistent with 
the ionization state found in Magellanic Stream sightlines 
with similar HI column densities (Fox et al. 2014). Fi¬ 
nally, deep 21-cm observations of the M81/M82 group reveal 
that relic gas from galaxy interactions can extend to cover 
large areas, possibly explaining the high Mgll covering frac¬ 
tions observed out to scales of < 200 kpc around luminous 
quasars. Tidal debris therefore represents a viable possibility 
for explaining observed Mgll absorption properties around 
quasars. We note, however, that it is not clear that tidal 
debris alone can explain the large fraction (30 — 40%) of 
Mgll absorption components found at radial velocities of 
|Au| > 300 kms“^. 


4.4 Future prospects 

The high covering fraction of Mg II absorption observed 
around quasars, correlation with luminosity, and complex 
kinematics observed in our SDSS sample are not consistent 
with the absorption expected from inactive galaxies with 
masses similar to typical quasar hosts. In addition, clus¬ 
tering measurements imply that high- and low-luminosity 
quasars are hosted by halos of similar mean masses (Shen 
et al. 2013b) implying that the correlation between cool gas 
covering fraction and quasar luminosity is not the result 
of simple mass scaling of the circum-galactic medium. To¬ 
gether, these observations imply that a substantial portion 
of the Mgll absorbing gas does not originate in the “nor¬ 
mal” halo gas of the quasar hosts. The high covering frac¬ 
tion, correlation with luminosity, and complex kinematics 
can be explained if the absorbing gas is the result of: (1) 
the GGM of neighboring galaxies on Mpc scales, (2) feed¬ 
back from luminous quasars, or (3) relics from the interac¬ 
tions thought to trigger luminous quasars. The first of these 


scenarios is in tension with the lack of correlation between 
quasar luminosity and clustering on Mpc scales (e.g. Berber 
et al. 2006; Shen et al. 2013b). The remaining two scenarios 
make distinct predictions that can be tested with additional 
observations. 

If a substantial fraction of the Mg II absorbing gas is the 
result of quasar feedback, then we expect the cool gas ab¬ 
sorption to be correlated with the presence of extended out¬ 
flows observed in [OIII] emission (e.g. Greene et al. 2012). 

If, on the other hand, the Mgll absorbing gas arises in 
relics from the mergers that are thought to trigger luminous 
quasars, then the absorption would be correlated with the 
presence of disturbed host morphologies and nearby tidal 
remnants. At z « 1, searching for such interaction signa¬ 
tures requires the high resolution imaging capabilities of the 
Hubble Space Telescope. 
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